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In this paper, two important factors which affect the pedestrian outflow at a bottleneck signifi- 
cantly are studied in detail to analyze the effect of an obstacle set up in front of an exit. One is a 
conflict at an exit when pedestrians evacuate from a room. We use floor field model for simulating 
such behavior, which is a well-studied pedestrian model using cellular automata. The conflicts have 
been taken into account by the friction parameter. However, the friction parameter so far is a con- 
stant and does not depend on the number of the pedestrians conflicting at the same time. Thus, we 
have improved the friction parameter by the frictional function, which is a function of the number 
of the pedestrians involved in the conflict. Second, we have newly introduced the cost of turning 
of pedestrians at the exit. Since pedestrians have inertia, their walking speeds decrease when they 
turn, and the pedestrian outflow decreases. 

The validity of the extended model, which includes the frictional function and the turning function, 
is verified by both a mean field theory and experiments. In our experiments, the pedestrian flow 
increases when we put an obstacle in front of an exit. The analytical results clearly explains the 
mechanism of the effect of the obstacle, i.e., the obstacle blocks pedestrians moving to the exit and 
decreases the average number of pedestrians involved in the conflict. We have also found that an 
obstacle works more effectively when we shift it from the center since pedestrians go through the 
exit with less turning. 

PACS numbers: 02.50.-r, 02.50.Ey, 34.10.+X, 45.70.Vn 



I. INTRODUCTION 

Pedestrian Dynamics has been studied vigorously in 
physics field over last decades [H, Si- Many 

microscopic 

models are developed such as the floor field (FF ) m odel 
@J JJ,Hi JMEII, the social force model 0,111 El, 
and the lattice-gas model d, EH EH , to simulate a crowd 
of pedestrians movement realistically. In addition to the 
simulations, there are also many pedestrian dynamics ex- 
periments [H, EE EH EH to study collective behaviors 
of pedestrians, while individual characteristics have been 
experimentally studied in psychology and physiology [2fj| . 

The behavior of pedestrians at an exit is focused on by 
many researchers in physics since it greatly affects the to- 
tal evacuation time in an emergency situation [ToL Il8l . 
Researchers of multi-robot systems have also studied such 
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evacuation [21| . In the Ref. 0, it is indicated that we 
can evacuate faster when an obstacle is put in front of 
an exit. In this paper, we have extended FF model for 
an analysis on the pedestrian outflow with an obstacle 
by introducing two important factors of pedestrian dy- 
namics at an exit . One is a conflict and the other is a 
turning. A conflict occurs when more than one pedes- 
trian move to the same place at the same time. When 
many conflicts occur, the pedestrian outflow, which is the 
number of pedestrians going through an exit with unit 
width per unit time, decreases, and the total evacuation 
time increases. Moreover, pedestrians must change their 
proceeding direction to the exit, when they go through 
it. Since pedestrians have inertia, although it is much 
smaller than that of vehicles, their walking speeds de- 
creases when they turn. This decrease in walking speed 
results in the decrease of the pedestrian outflow. 

In the FF model, which is a pedestrian dynamics model 
using cellular automata, conflicts are take into account by 
friction parameter Q , which describes clogging and stick- 
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FIG. 1: (a) A schematic view of an evacuation simulation by the FF model. The person-shaped silhouettes represent 
pedestrians; the letters E and O represent the exit cell and obstacle cells, respectively. Pedestrians proceed to the exit by one 
cell at most by one time step, (b) Target cells for a pedestrian at the next time step. The motion is restricted to the Neumann 
neighborhood in this model. 



ing effects between pedestrians. Since the friction param- 
eter is a constant parameter, the strength of clogging and 
sticking does not depend on the number of pedestrians 
involved in the conflict, fn reality, however, it is more 
difficult to avoid a conflict when three pedestrians move 
to the same place at the same time than two pedestrians 
move to the same place at the same time. The inertia 
of pedestrians is considered in Ref. @, HI, HH , however, 
this inertia effect is the suppression of quick changes of 
the direction of the motion, i.e., pedestrians try to keep 
walking in the same direction. The inertia effect, which 
should be introduced at an exit, is not the suppression of 
quick changes, but the decrease in walking speed, which 
has not been considered in the previous models. 

Therefore, we newly introduced the frictional function 
and the turning function. Frictional function is a func- 
tion of the number of pedestrians involved in the conflict, 
which reproduces the difference of the strength of clog- 
ging against it. Turning Function is a function of the 
angle of deviation from the former direction of motion, 
which estimates the decrease in walking speed quantita- 
tively by turning. These functions make it possible to 
describe the pedestrian behavior around an exit more 
precisely and realistically. The mechanism of the effect 
of the obstacle is also explained successfully by the fric- 
tional and turning functions. 

We have also performed the experiments of evacuations 
and verified that the outflow obtained by the extended 
model agrees with the experimental data well. The out- 
flow increases when an obstacle is put at the proper po- 
sition in front of the exit as we have expected from our 
theoretical calculation. 

We would like to emphasize that our study is based on 
firm theoretical analysis. Vehicle traffic has been theoret- 
ically studied by extending asymmetric simple exclusion 
process [13, HH, HH, however, pedestrians dynamics 
has been hardly done because of the complexity of the 
rules of motion and the two-dimensionality. Therefore, 
theoretical analysis on pedestrian dynamics is also a new 



challenge. 

This paper is organized as follows. In Sec. [D] we briefly 
review the FF model. The two new functions, i.e., the 
frictional function and the turning function, are intro- 
duced in Sec. |TTT] and Sec. IIV1 respectively. In Sec. |V] 
several kinds of lattice for FF model is considered. We 
obtain the analytical expression of the average pedestrian 
outflow by using the cluster approximation in Sec. IVII 
In Sec. IVIII and Sec. IVIII1 we show the result of the 
evacuation experiments in lines and the evacuation ex- 
periment with an obstacle, respectively, and investigate 
the mechanism of increasing outflow by obstacle. Section 
IIXI is devoted to summary and discussion. 



II. FLOOR FIELD MODEL 

A. Floor Field 

We consider a situation that every pedestrian in a room 
moves to the same exit. The room is divided into cells 
as given in Fig. [T] (a). Man shaped silhouettes repre- 
sent pedestrians, an alphabet E and alphabets O repre- 
sent the exit cell and obstacle cells, respectively. Each 
cell contains only a single pedestrian at most. Every 
time step pedestrians choose which cell to move from 5 
cells for simplicity: a cell which the pedestrian stands 
now = (0,0)] and the Neumann neighboring cells 

= (0,1), (0,-1), (1,0), (-1,0)] (Fig. [H (b)). Two 
kinds of FFs determine the probability of which direc- 
tion to move, i.e., Static FF Sy, which is the shortest 
distance to the exit cell (Fig. [5]), and Dynamic FF Dij, 
which is a number of footprints left by the pedestrians 
4]. Dynamic FF plays role in mimicking the long-raged 
interaction among pedestrians to short-ranged one. The 
transition probability pij for a move to a neighboring cell 
(i,j) is determined by the following expression, 

Pij = N&j exp(-k s Sij + kdDij). (1) 
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FIG. 2: Static floor field constructed by the exit E. The 
numbers in each cell represent the Euclidean distances from 
the exit cell, which are calculated by Dijkstra method [28j ]. 



k = 2 



l-ju 



14 



or 



FIG. 3: The way of solving conflicts. In a conflict situation 
(k > 2), movement of all involved pedestrians is denied with 
probability fi. One of them is randomly allowed to move to 
the desired cell with probability 1 — fi. 



Here the values of the FFs Sij and at each cell (i,j) 
are weighted by two sensitivity parameters k s and ka with 
the normalization N . There is a minus sign before k s 
since pedestrian move to a cell which Static FF decreases. 
£ij returns for an obstacle or a wall cell and returns 1 
for other kinds of cells. Note that in our paper a cell 
occupied by a pedestrian is not regarded as an obstacle 
cell, thus it affects the normalization N. In the following, 
we ignore the effect of Dynamic FF, i.e., k& — 0, since we 
study the pedestrian behavior at an exit, where there is 
only short-ranged interaction. 



B. Bottleneck Parameter 

We use bottleneck parameter (3 6 [0, 1], which was in- 
troduced in Ref. [l0(, when we consider the congested 
evacuation situation at an exit. Then, the transition 
probability of pedestrians who occupy one of the neigh- 
boring cells of the exit cell is transformed from pij to 
(j>ij)p as follows: 



(Pij)p = 
(Po,o)/3 



PPi. 



[(i,j)^(Q,0)], 



Ppo,0 + (1 - Ei i PPij) = PPo,o + (1 - P) 



(2) 



fi controls the velocity of the pedestrians who are at the 
neighboring cells of the exit. When k s is large, the tran- 
sition probabilities of pedestrians at neighboring cells of 
the exit are approximated by using (3 [10|]. This sim- 
plification enables us to analyze the pedestrian behavior 
theoretically. 



C. Friction Parameter 

Due to the use of parallel dynamics it happens that 
two or more pedestrians choose the same target cell in 
the update procedure. Such situations are called conflicts 
in this paper. To describe the dynamics of a conflict in a 



quantitative way, friction parameter \i was introduced in 
Refs. This parameter describes clogging and sticking 
effects between the pedestrians. When we denote the 
number of pedestrians moving to the same cell at the 
same time as k £ N, the friction parameter is not used 
in k = 1, but applied in k > 2 as follows: 







(fc = i), 

(fc>2), 



(3) 



where M is the probability that all the pedestrians in- 
volved in the conflict remain at their cell. <pfj,(l) — since 
there is no conflict in k = 1. In k > 2, the movement 
of all involved pedestrians is denied with probability /i, 
therefore, the conflict is solved with probability 1 — fi, 
and one of the pedestrians is allowed to move to the de- 
sired cell (Fig. [31). The pedestrian which actually moves 
is then chosen randomly with equal probability. In a sit- 
uation with large [i, pedestrians are competitive and do 
not give way to others. Thus they hardly move due to 
the conflict between them. Contrary in a situation with 
small /i they give way and cooperate each other. 



D. Update rules 

The FF model which ignores the effect of the Dynamic 
FF consists of the following two steps per unit time step. 

1. Calculate each pedestrian's transition probability 
by (P) and ^ . Transition probabilities of pedestri- 
ans who stand on the exit cells, i.e., the probabili- 
ties of getting out from the room, is a € [0, 1]. 

2. Move pedestrians based on the calculated transi- 
tion probability with parallel update. If there are 
cells which are possibly occupied by more than one 
pedestrian, solve conflicts by the means of Sec. Ill CI 
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FIG. 4: The number of conflicts in a competitive situation, 
i.e., k 3 — 20, /3 = 1, and u = 0.6. We simulated in the 11x11 
cells' room with the 1 cell's exit for 101,000 time steps and 
accumulated the value of 1,001 to 101,000 time steps. The 
exit cell is at (6, 10) and Neumann Neighborhood is used. 



FIG. 5: Values of the frictional parameter and the frictional 
function against k, which is a number of pedestrians who move 
to the same cell at the same. We see that <^ is constant in 
k > 2, while 4>i increases gradually and saturate in 1. (p^ 
and 4>( become large according to the increase in fi and 
respectively. 



III. FRICTIONAL FUNCTION 

Figure 0] shows the number of conflicts, which is a re- 
sult of a evacuation simulation using Neumann Neigh- 
borhood in 100,000 time steps. In this simulation pedes- 
trians enter into the room with probability 1 from the 
boundary cells of three sides of the square room that do 
not include the exit cell, therefore the room is always 
filled with pedestrians [10| . At the exit cell, there are 
69,385 conflicts in 100,000 time steps, i.e., the probabil- 
ity of a conflict at there is approximately, 70%. Thus, 
the conflict is the important factor when we consider the 
evacuation from the narrow exit, and the friction param- 
eter has played a role to represent it. However, it is 
constant in k > 2 and does not take account of the dif- 
ference of the number of the pedestrians involved in the 
conflict. In reality, an impact of a conflict is stronger 
when three pedestrians conflict each other than when 
two pedestrians conflict each other. In the simulation, 
85% of the conflicts occur in the room are the conflicts 
between two pedestrians, and 15% of them are the con- 
flicts among three pedestrians. When we consider the 
conflicts at the exit cell, 34% are between two pedestri- 
ans and 66% are among three pedestrians. Thus, it is 
impossible to approximate as "Number of the conflicts 
between two pedestrians" <C "Number of the conflicts 
between three pedestrians" or the opposite way, and we 
cannot consider the only one kind of conflict. 

Therefore, we newly introduce the frictional function 
4>q G [0, 1], which is a function of the number of pedestri- 
ans moving to the same cell at the same time k in general 
as follows: 

Mk) = 1 - (1 - () k - k((l - C)^ 1 (*>1). (4) 
4>r represents the probability of the unsolved conflict and 



is defined by considering the psychological effect of pedes- 
trians. £ £ [0, 1] is an aggressive parameter, which is a 
probability of not giving way to others when more than 
one pedestrian move to the same cell at the same time. 
The term (1 — £) fe in the expression ([4} is the probability 
that all pedestrians involved in the conflict try to give 
way to others. The term fc£(l — £) 1 is the probabil- 
ity that only one pedestrian does not give way to others 
while the others do. By subtracting the two terms, which 
are the probabilities of resolving the conflict, from 1, we 
obtain the frictional function <f>^. Therefore includes 
the effect of motivations of pedestrians. 

We compare the characteristics of cf>^, i.e., frictional 
function and which we call frictional parameter in 
the following. <f>^ satisfies the two conditions: 



0C(1) = O, 

0c(°°) = 1 



(C^o). 



(5) 



The former equation means that there is no conflict when 
only one pedestrian move to the target cell. The latter 
equation means that no one can move to the same cell 
when greatly many pedestrians move to there at the same 
time since pedestrians have some finite volume, (j)^ sat- 
isfies the former equation, however if /i ^ 1, it does not 
satisfy the latter one, which is a natural condition near 
the exit. 

Figure [5] is the plots of </) M and <j>^ as a function of k. 
4>^ is constant in k > 2, while 4>( gradually increases as 
k increases and reflects the difference of the strength of 
clogging against k. We compare the validity of 4>^ and 
<f>Q by applying them to the experimental results in the 
Sec. IVTfland lVlTIl 

The advantage of the frictional function is that its 
value depends on the number of pedestrians involved in 
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FIG. 6: Values of the turning function r against 9, which is 
an angle of deviation from the former direction of motion. It 
decreases as 8 increases and its slope become sharp when r\ 
increases. 



a conflict, but does not depend on the type of the lattice, 
which we adopt. Thus, we can use them in every discrete 
model of self-driven particles with parallel update. We 
would also like to emphasize that the frictional function 
is not only effective for considering a narrow exit, i.e., an 
exit with one cell, but also a wide exit, i.e., an exit with 
plural cells, since the friction parameter is also ap plie d 
to calculate an outflow through a wide exit in Ref. [ig |. 



IV. TURNING FUNCTION 

Since pedestrians have inertia, their walking speeds 
decrease when they turn and change their headings. 
While the inertia effect of vehicles is considered in Ref. 
pEL [20! . |3C| as a slow-start rule, that of pedestrian dy- 
namic is introduced to the FF model in Ref. 0, [HI, HH • 
However, these inertia effects are not appropriate for es- 
timating the decrease in walking speed. 

The inertia parameter introduced in Ref. Q is too 
simple for analyzing the pedestrian behavior at the exit 
in detail. When the strength of inertia increases in the 
model, pedestrians try to move the same direction as they 
moved in the previous step. Therefore, the pedestrian 
move to the exit from the lateral direction often cannot 
stop at the exit and pass in front of it. Moreover, the 
strength of inertia does not depend on the angle between 
the direction of the movement in previous step and that 
in present step. For instance, the probability of turning 
90° and that of turning 180° are same. In reality, how- 
ever, it becomes more difficult to turn as the angle of 
deviation from the former direction becomes large. It is 
also verified by the experiments in Ref. [2(| • 

In the F. A. S. T. model [H, [H[, an inertia effect is 
studied in detail. The strength of the inertia effect is 
a function of the angle of deviation from the former di- 
rection and calculated by considering centrifugal forces. 
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FIG. 7: A schematic view of the turning pedestrian. The 
pedestrian, who is depicted with circled-arrow in the figure, 
moved to the cell (0,0) from (0, —1) in the time step t — 1. 
If he/she try to move to the cell (1,0), he/she has to turn 
0i,o(= —90°) to change his/her heading, and if he/she try to 
move to the cell ( — 1, 1), he/she has to turn 0_i,i(= 45°) to 
change his/her heading. Oij changes its value according to the 
former direction of the movement. 



However, this inertia effect is not used to decide the walk- 
ing speed but only the moving direction. 

Thus, we introduce the turning function t, which rep- 
resents the effect of decrease in walking speed by turning 
as follows: 



r(0)=exp(-J7|0|) 



(6) 



where 8 G [tt, —tt] is the angle of deviation from the for- 
mer direction of motion, and 77 > are the inertia co- 
efficient in turning, which represent the strength of the 
inertia. We decide to use an exponential function for the 
turning function since it is simple and also matches to 
the transition probability ((T|) in the FF model. An expo- 
nential function is also easy to analyze theoretically. We 
see that r(0) = 1 for arbitrary 77 from ^ since there is no 
decrease in walking speed when pedestrians do not turn. 
Figure [5] shows that the values of the turning function, 
which decrease almost linearly according to the increase 
of the angle of deviation 9 since rj is not large. We verify 
that 77 for real pedestrians is as small as 0.1 in Sec. IVIII 
and I VIIH The turning function enables us to calculate 
more realistic outflow. 

(Pij)r, which is the transition probability including the 
effect of turning by r, is described as, 



(Pij)r = T{9ij)pi 



Km) ^(o,o)], 



(po,o)t = t(0)po,o + 1 - ^2 T (8ij)Pi. 



(7) 



where pij is given by Jl}, an d 0-y is the angle between 
the former direction of the motion and the direction of 
the cell as in Fig. [7] 
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FIG. 8: Schematic views of the pedestrian cluster at the exit cell for different neighborhood, (a), (b), and (c) are the case 
of Neumann neighborhood, Triangle-lattice neighborhood, and Moore neighborhood, respectively. Circled-arrows represent 
pedestrians and their facing direction. In the first row, target cells for a pedestrian are depicted. In the case of Neumann 
neighborhood, pedestrians can move only vertical and horizontal directions, however, they can move to diagonal ways in the 
case of Moore neighborhood. Pedestrians can move to six directions in the case of Triangle-lattice neighborhood; however 
they cannot go straight in the vertical way. In the second row, the schematic views of pedestrian cluster at the exit cell and 
its neighboring cells are depicted. The incident angle 8 e is defined as an angle between the line l C enter and direction of the 
pedestrian movement as in the figure (a). For instance, pedestrian-3 tries to move to the exit cell from the right neighboring 
cell of the exit cell in the figure (a). Thus his/her angle of incidence 8 e equals to —90°. n e , i.e., the numbers of the neighboring 
cells of the exit cell, and ranges of fe e , i.e., the number of pedestrians moving to the exit at the same time are also described in 
the figure. The minimum and the maximum value of k e are and n e , respectively. 



V. NUMBER OF THE NEIGHBORING CELLS 
AND THE INCIDENT ANGLE 

In this section, we consider the number of the neighbor- 
ing cells of the exit cell and the incident angle to the exit. 
We denote them as n e and 9 e in the following. Since FF 
model adopts Neumann Neighborhood, there are three 
neighboring cells of the one-cell exit, i.e., n e = 3 (Fig. [5] 
(a)), thus, fc e , which is the number of pedestrians moving 
to the exit at the same time, is between and n e {= 3) 
and, there is possibility of occurring a conflict between 
two pedestrians and three pedestrians at the exit. We 
suppose that all pedestrians at the neighboring cells face 
to the exit cell in the evacuation. Then, the incident 
angle 8 e is defined as an angle between the vertical line 
perpendicular to the wall (l cen ter in Fig- [HI (a)) and the 
direction of the pedestrian movement. Therefore, the in- 
cident angle of the pedestrian at the cell 2 in Fig. [5] (a) 



is 0°, and those of pedestrians at the cell 1 and 3 are 90° 
and —90°, respectively. They have to turn 9 e at the exit 
cell to change their orientation and go through the exit 
(Fig. 0. 

Here, we consider two other neighborhoods in the FF 
model. One is the Triangle-lattice neighborhood, which 
is obtained from setting up cells on the node of the tri- 
angle lattice. In the case of the Triangle-lattice neigh- 
borhood, there are four neighboring cells of the exit cell, 
i.e., n e = 4, thus a conflict between two, three, and four 
pedestrians occur. The incident angles are calculated as 
in Fig. [5] (b) since we see equilateral triangles by con- 
necting the center of each cell. Note that the Fig. [H] (b) 
is just a schematic view, and the each cell is depicted 
as approximated size. The situation as in Fig. [8] (b) is 
realistic for the congested evacuation since pedestrians 
often stand between the two former pedestrians to see 
their way clear between the two heads. The Triangle- 
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FIG. 9: Schematic view of the pedestrian-1 going through the exit, (a)— »(b): Move to the exit cell, (b)— >(c): Turn 90° at the 
exit cell, (c)— »(d): Get out from the room. The transition probability from (a) to (b) is represented by p(l — » E) l[8]). and that 
from (b) to (d) is described by p(E — > out) ([9]). Hence, p(E — > out) includes both turning and getting out. 



lattice neighborhood also reproduces the hexagonal close 
packing of pedestrians as the two dimensional granular 
flow. The other is the Moore neighborhood as in Fig. 
IS] (c). Pedestrians have eight choices when they move 
and there are five neighboring cells of the exit cell, i.e., 
n e = 5, thus there is a possibility of occurring a conflict 
between two to five pedestrians at the exit. The incident 
angles are as in Fig. [5] (c) . The difference between the 
Neumann neighborhood and the Moore neighborhood in 
the FF model is studied in Ref. [3lj |. 

When we consider the pedestrian dynamics in the low 
density case, where there are a few pedestrians, the num- 
ber of neighboring cells does not greatly affect on the 
number of conflicts, however, it becomes significantly dif- 
ferent according to the neighborhood type in the high 
density case. Since we study an evacuation situation, the 
density at the exit is always high, thus we verify which 
neighborhood is the most suitable for analyzing evacua- 
tion by our experiments in Sec. I VIII 



VI. ANALYTICAL EXPRESSION FOR THE 
AVERAGE PEDESTRIAN OUTFLOW USING 
THE CLUSTER APPROXIMATION 

A. Approximate Transition Probability 

We study the theoretical expression for the average 
pedestrian outflow through an exit in an extremely con- 
gested situation for arbitrary kinds of neighborhoods by 
cluster approximation in the FF model. We consider a 
room with an exit with one-cell, and focus on the exit cell 
and its neighboring cells 1 to n e in Fig. [TO] The num- 
ber of the neighboring cells of the exit is described as 
arbitrary n e 6 N to obtain the expression of the outflow 
for arbitrary kinds of neighborhoods. The case n e = 3, 
n e = 4, and n e = 5 correspond to the case Neumann 
neighborhood, Triangle lattice neighborhood, and Moore 
neighborhood, respectively, as in Fig. [5] 

We assume k s — > oo since pedestrians near the exit try 



to evacuate in the shortest way. Then, a pedestrian at 
neighboring cell m 6 [l,n e ], which we call pedestrian- m 
in the following, has only two choices, i.e., moving to the 
exit cell with the probability p(m — > E) or stay at his/her 
cell p(m — > to), which are described as 

| p(to-s-E) = (py)/s|(*,j)=E , g . 
\ p(m ->• to) = (po,o)/3 (k s -> oo), 

respectively. 

A pedestrian who is at the exit cell also has only two 
choices, i.e., getting out with the probability p(E — * out) 
or stay at the exit cell with the probability p(E — > E), 
which are described as follows: 

f p(E -> out) = ar{8 e . m ), 
\p(E-fE) =l-ar(6 e , m ), 

where a is the probability of getting out from the room, 
which does not include the effect of turning (Sec. Ill Dp . 
and B e . m is the angle of incidence when pedestrians move 
from the cell m . We denote 6^ m as m in the following 
for simplicity as in Fig. [TO] 

p(E — > out) represents both turning at the exit cell 
and getting out from the room (Fig [9]). The turning 
function ^ is introduced to transition probability of the 
pedestrians at the exit cell since they need to change their 
headings quickly to face toward the exit (Fig. [S](b)— >(c)), 
while we neglect the effect of turning at the neighboring 
cells of the exit for simplicity since we assume that all 
pedestrians at the neighboring cells are directed to the 
exit cell in the evacuation situation. 



B. Mathematical Formulation of 
the Pedestrian Outflow 

Each neighboring cells of the exit has two states, i.e., 
occupied by a pedestrian and vacant, and the exit cell 
has n e + 1 states since we consider the headings of the 
pedestrians to introduce the turning effect. Therefore, 




n e pedestrians at n e neighboring cells 

FIG. 10: Schematic view of the cluster at the exit cell and 
its neighboring cells. We consider the congested evacuation 
situation, thus all n e neighboring cells are occupied by pedes- 
trians. They try to move to the exit cell and conflict each 
other. (3 and ar(9 ra ) represent the transition probability to 
the exit cell (JSj , and the transition probability getting out 
from the exit cell (SJ, respectively. 6 m (m £ [l,n e ]) is the 
incident angle when the pedestrian-m moves to the exit cell. 



there are (n e + l)2 ne states when we focus on the cluster 
at the exit cell and its neighboring cells. We assume 
that pedestrians come into the neighboring cells of the 
exit with probability 1 since we consider the congested 
evacuation. Then, the number of the states of the cluster 
reduces to 2n e + 1. Figure [Til shows the reduced states of 
the cluster in the Neumann neighborhood case (n e = 3). 
There is a state that the exit cell is vacant, i.e., (SO), 
and there are two states that the exit cell is occupied 
by pedestrian-m, i.e., (Sm-A) and (Sm-B), for each m(= 
1,2,3) in the Fig. [JU 

We need the probability distribution of the exit-cell 
state to calculate the pedestrian outflow through the exit; 
however, do not need that of neighboring cells. Thus, we 
combine the state (Sm-A) with (Sm-B) and consider the 
state transition among n e + 1 exit-cell states. Figure 
[T2l is the state transition diagram of the exit cell in the 
Neumann neighborhood case (n e = 3). We find that the 
diagram is symmetrical with respect to the (SO). ar(0 m ) 
is the probability of getting out from the exit cell given 
by ©. r(n e ) represents the probability of one pedestrian 
move to the exit cell by solving a conflict when there are 
n e neighboring cells of the exit. As we describe in Sec. 
Ill C[ the pedestrian which actually moves is chosen ran- 
domly with equal probability when the conflict is solved. 
Thus, the transition probability of (SO) to (Sm) is given 
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FIG. 11: Schematic view of 7(= 1 + 2x3) states of the cluster 
in the Neumann neighborhood case (n e = 3). (SO) represents 
the state that the exit cell is vacant, and (Sm) represents 
the state that the exit cell is occupied by pedestrian-m (m = 
1,2,3). The pedestrian-1, 2, and 3 are circled — >, circled j, 
and circled in the figure, respectively. The States (Sm-A) 
and (Sm-B) in each (Sm) are combined when we calculate the 
pedestrian outflow. 



by r(n e )/n e . An explicit form of the r(n e ) is obtained in 
the following. 

Now, we start to calculate the average pedestrian out- 
flow. First, the probability of k e £ N pedestrians trying 
to move to the exit cell at the same time is described as 



b(k e ) 



(10) 



by using ([5]). Next, r(n e ), which is the probability of one 
pedestrian succeeds to move to the exit cell, is obtained 
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FIG. 12: State transition diagram of the cluster in the Neu- 
mann neighborhood case (n e — 3). r(n e ) is the probability 
that one of the pedestrian at the neighboring cells succeeds to 
move to the exit cell by solving conflicts (111[) . and ar(6 | m ) is 
the probability of getting out from the exit cell © . Since the 
diagram is symmetrical with respect to the (SO), we success- 
fully calculate the average pedestrian outflow with arbitrary 
number of the neighboring cells n e . 



we obtain the stationary (t — > oc) solution 



7Too(0) 



7roo(m) 



r(n e ) 



i 

y -4— 

n e a f— ' T(drh) 



(14) 



r(n e ) 



-ttoo(O) (m G [1, n e ]). 



n e aT(9 m ) 

Thus, the number of pedestrians who can evacuate from 
an exit with one cell's width per one time step, i.e., the 
average pedestrian outflow through the exit, is described 
as 

n e 

(qtheo(n e ,6i,...,6 ne )) = ^ aT(9 m )n ao (m) 



^(«e)7Too(0) 



-i -1 



1 | _L_ 1 

r(n e ) n e a f-' r(^) 

m— 1 

(15) 



where (x) represents the sample average of x. 



C. Four Formulations of the Pedestrian Outflow 



as 



in e )= £ [{l-4>(k e )}b(k e )}, 



(11) 



where cj) is the frictional parameter or function given by 
© or (U]), respectively. We define 7r t (0) as the probability 
that the cluster is at the state (SO), i.e., a pedestrian is 
not at the exit cell at time step t and 7r t (m) (m € [1, n e \) 
as the probability that the cluster is at the state (Sra), 
i.e., a pedestrian who was at the neighboring cell m is at 
the exit cell. Then, the master equation of the cluster is 
described as follows: 



Wo) 



7r t+ i(n e ) 

1 - r(n e ) ar(6*i) 
r(n e )/n e 1 — ar(di) 



"(n e )/n e 







ar(9 n J 



1 - aT(9 ne 



Vt(0)" 

Tt(l) 
TTt(^e) 



(12) 



where ar(9 e ) represents the probability of getting out 
from the exit ©. By using (|12[) with the normalization 
condition 



(13) 



There are two choices to represent conflicts, i.e., using 
frictional parameter </> M or frictional function 4>^. We also 
consider whether the turning effect is included or not 
for each outflow. Therefore, we obtain four kinds of the 
theoretical pedestrian outflows as follows: 

• (<?a»o) = (9thco)U=0„,» ) =O 

The simplest formulation of the pedestrian outflow, 
which includes the effect of conflict by the frictional 
parameter and neglects the turning effect. The re- 
sults obtained in Ref. [fc| is recovered by substi- 
tuting n e = 1,2,3 as follows: 



(9moW) 

<?M)(2)) = 

(«mo(3)> = 
1 



a/3 
CV + /3' 



a + 2/3 - (1 + /i)/3 2 



a 



a + 3f3 - 3(1 + n)(3 2 + (1 + 2 A i)/3 3 



(16) 



(17) 



(18) 



• (iCo) = (<Ztheo)U=^,»7=0 

The pedestrian outflow which includes the effect 
of conflicts by the frictional function and does not 
consider the effect of turning. This formulation is 
obtained in Ref. 1321 as 



(%o(n e )) 



1 1^1 

-r-; + —/ ~ 

r(n e ) n e ' a 

. m— 1 

ar(n e ) 
a + r(n e ) 



(19) 
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FIG. 13: (a) Average pedestrian outflows (q^o) and (q^o) 
as functions of the number of the neighboring cells of the 
exit n e in the case (5 = 1, a = 1. (q^o) does not change in 
n e > 2, while (q^o) keeps decreasing in the all region in n e > 1. 
(b) Average pedestrian outflows (500) and (qov) f° r different 
evacuation types in the case /3 = 1, a = 1. The explanation 
of the three evacuation types are described in Fig. [T3] We see 
that (goo) does not change according to the three evacuation 
types, whereas (qov) does change. The slope becomes sharp 
when 77 increases. 



(<W - (9theo) U=<^„ 

In this pedestrian outflow formulation, conflicts 
and turning are considered by the frictional param- 
eter and the turning function, respectively. 



(SCn) = (<7thoo)U=0 c 

The pedestrian outflow, which includes the effect 
of conflicts and turning by the two new functions 
introduced in this paper, i.e., the frictional func- 
tion (j)^ and the turning function r. 



In the Sec. IVIII and lVIIIl we verify which formulation is 
the most suitable for calculating the realistic pedestrian 
outflow. 
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FIG. 14: Schematic views of the three evacuation types in 
the case (n e = 3). The neighboring cells are denoted as 1, 2, 
and 3. The figures described below the schematic views are 
the incident angles (#1, 62, 63). The names of the evacuation 
types are decided from the directions from which pedestrians 
come into the exit cell. 



D. Comparison of the Theoretical Outflow 



Figure [13] (a) shows the average pedestrian outflow 
(<Z/jo) an d (9^0)1 which are outflows that includes the ef- 
fect of conflicts but does not include the effect of turning. 
We see the difference between the effect of the frictional 
parameter M and the frictional function ip^ from the fig- 
ure. When we use frictional parameter, the outflow does 
not change as n e , which is the number of neighboring cells 
of the exit cell, increases in n e > 2, while it continues to 
decrease when we adopt the frictional function. Since it 
is natural to consider that the strength of the clogging 
and sticking become stronger according to the increase of 
the number of pedestrians involved in the conflicts, (<?£o) 
seems more realistic than (g M o)- We obtain the decrease 
of the outflow in n e > 2 in the urgent evacuation situa- 
tion (f3 = 1) by the frictional function for the first time 
in this paper. 

Figure [13] (b) shows the average pedestrian outflow 
(qoo) and {qo v )- In these cases, the effect of conflicts is 
neglected by substituting for fi or £, and we see how the 
effect of turning influences on the outflow. The schematic 
views of the three evacuation types, which are V-type, L- 
type, and T-type are depicted in Fig. Q3J n e — 3 in 
all three types, however, the incident angles of the three 
pedestrians 9\, 8 2 , and 9 3 are different. When we neglect 
the turning effect ((?oo))i there is no difference of the out- 
flow among the three evacuation types. However, when 
we consider the effect of turning by the turning func- 
tion ((qori)), the outflow decreases as the incident angles 
increase. The difference of the pedestrian outflow accord- 
ing to the difference of the way of pedestrians gathering 
around the exit is observed since we have introduced the 
significant factor in the evacuation through the narrow 
exit, i.e., the decrease in walking speed by turning, by 
the turning function. 
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FIG. 15: (a) Schematic view of the room used in the experiments. The width of the door is 50 cm with a locker at both 
sides, (b) Schematic views of nine conditions of the experiments. In the case from (A) to (G), pedestrians evacuate in lines. 
They are forbidden to pass the former pedestrian. The figures described by the silhouettes of pedestrians are the approximated 
incident angles used in theoretical calculation. The participants of the experiments are not exactly go to the exit along the 
angle of incidence described in the figures. The way of the evacuation is same in (H) and (I), however, the initial conditions 
are different as in the figures. 



VII. EXPERIMENT AND ANALYSIS 1 
EVACUATION IN LINES 

A. Experiment 

We have performed the evacuation experiments to ver- 
ify the relation between the two new factors and the aver- 
age pedestrian outflow. The schematic view of the room 
is described as in Fig. [15] (a). The width of the exit 
is 50 cm. There are eighteen participants for the experi- 
ment, who are all men. The experiments started when we 
clapped our hands and finished when the all pedestrians 
evacuated from the room. Nine kinds of conditions are 
put into practice (Fig. \T5\ (b)). In the case (H) and (I) 
pedestrians could move as they want after the evacuation 
started, however, in the other cases pedestrians had to 
follow the former pedestrian, i.e., they were prohibited 
from putting the queues into disorder. Each experiment 
was performed two or three times, which is sufficient for 
pedestrian dynamics experiments. When we perform the 
same experiments many times, pedestrians memorize and 
study the situation, which greatly affects the experimen- 



tal results. 

The outflows are calculated as follows: 



(<?exp) 



w(tj - U) ' 



(20) 



where i and j are the orders of the first and the last 
pedestrian used to calculate the pedestrian outflows, and 
ti and tj are the evacuation times of the i th pedestrian 
and the j th one, respectively, w represents the width of 
the exit. 

We decide i — 4, j = 18 for the case (A) to (D) and 
i = 1, j = 15 for the case (E) to (H). In the former 
cases, pedestrians move very fast at the starting time 
since there is no one ahead of the first pedestrian and the 
density is very low. Thus, we need to wait for the stable 
state. By contrast, in the latter cases, the evacuation 
becomes saturated in the stable state immediately since 
there are interactions among the pedestrians from the 
start. However, in the end of the evacuation, the number 
of pedestrians decreases and the interaction among more 
than two pedestrians is hardly observed. Consequently, 
the situations become totally different from the stable 
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TABLE I: The pedestrian outflows of the experiments. 
"Case" in the table is corresponding to the case in Fig. 1151 (b). 
n e is the number pedestrian lines in the experiments, which 
represents the number of the neighboring cells in the theo- 
retical calculation. e — (9i, . . . ,0n B ) is the approximated 
incident angles of the pedestrians. The pedestrian outflows 
(g C x P ) [persons/(m- sec)] in the table are the average of the N 
experiments, where N is the number of the experiments de- 
scribed in the fifth column of the table, i and j are the orders 
of the first and last pedestrians which used to calculate the 
experimental outflow by (|20[) . 



Case 


n e 


e 


(<7cxp) 


N 


i 


j 


(A) 


1 


(0°) 


2.62 


3 


4 


18 


(B) 


2 


(30°, 30°) 


2.81 


3 


4 


18 


(C) 


2 


(0°,90°) 


2.69 


3 


4 


18 


(D) 


2 


(90°, 90°) 


2.59 


3 


4 


18 


(E) 


3 


(45°,0°,45°) 


2.69 


2 


1 


15 


(F) 


3 


(90°,0°,90°) 


2.59 


3 


1 


15 


(G) 


4 


(90°, 30°, 30°, 90°) 


2.51 


2 


1 


15 


(H) 






2.51 


3 


1 


15 


(I) 






2.53 


2 


4 


15 



states. Therefore, we use the evacuation time of first and 
fifteenth pedestrian for the calculation. In the case (I), it 
takes some time to achieve the stable state since pedes- 
trians are not gathering around the exit at the beginning 
of the evacuation. In the end of the evacuation, the same 
phenomenon as the experiment (E) to (H) is observed. 
Thus, we determine i — 4 and j — 15 for the case (I). 

The results of the experiments are described in Tab. 
HJ There are four remarkable points in the table. First, 
we see that the pedestrian outflow decreases in n e > 2 
by comparing the case (B), (E), and (G). This indicates 
that when the number of the pedestrian-lines increases, 
the outflow decreases since the possibility of the conflicts 
increases as in Fig. [SJ Second, we observe that the out- 
flow decreases with the increase of the components of 9 e 
when n e is constant by comparing (B), (C), and (D) or 
(E) and (F). We have verified that turning at the exit 
does decrease the outflow. Third, the maximum outflow 
is attained at n e = 2. See the Fig. [16j which is a snap- 
shot of the experiment (B). We see that pedestrians in the 
left line and the right line go through the exit one after 
the other. This phenomenon is called zipper effect 
when it occurs, the outflow increases since the number of 
conflicts between pedestrians decrease dramatically. In 
the traffic flow, the alternative configuration, which leads 
to the smooth merging as the pedestrian evacuation in 
the two-line case, is achieved by the compartment line 
[HI]. By contrast, pedestrians are intelligent enough to 
avoid conflicts by their selves. They do not need any 
help to achieve the alternative configuration. Forth, the 
pedestrian outflow in the case (G) (n e = 4), (H), and (I) 
are similarly smaller than the other cases. This indicates 
that there are approximately four pedestrians at the exit 
in the normal evacuation situation when the width of 



Exit 

50 cm 




FIG. 16: Snapshot of the experiment (B) (n e = 2). The 
pedestrians in the left line and the right line go through the 
exit one after the other without any instructions. 



the exit is 50 cm. Consequently, we have found that the 
Triangle-lattice neighborhood is the most suitable for the 
normal evacuation. 



B. Theoretical Analysis 

We calculate the average errors between the theoretical 
outflows and the experimental outflows as follows: 



ZXv2thco) ~ (9cxp)) 2 



(Number of the kinds of the experiments) 



(21) 



and verify the effect of the frictional function and the 
turning function. First, we define the length of the side 
of square cell Al and the time of one time step At as 



Al = 0.5 [m/cell], 
At = 0.3 [sec/step]. 



(22) 



Second, we determine a and (3. In the situation (A) to 
(I), the pedestrians do not change their walking speeds 
instantaneously after they succeed to move to the exit 
cell. Therefore, we assume a = (3. Then, we calculate 
from the result of the experiment (A) as follows: 



(<ZthcoK =1,01 = 0°)) = 2.62 x Al x At, 
(3 = 0.79. 



(23) 



Note that when n e = 1, the frictional parameter, func- 
tion, and turning function do not affect on the pedestrian 
outflow, so that we obtain common (3 for the four kinds 
of outflows. 

Finally, we decide the parameters /i or £, and r/ for 
four theoretical average pedestrian outflow introduced in 
Sec. I VI CI from the method of least squares. The number 
of the pedestrian lines is assumed as n e = 4 in the case 
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TABLE II: Parameters and the average errors for the experi- 
ment (A) to (I) from the method of least squares. We see that 
the average errors calculated by (|21|l are large in three cases, 
i.e., (qfio), (<?co}> an d (<7mij)- However, if the both frictional 
and turning functions are applied to the outflow, i.e., (q( v ), 
the error becomes smaller than half of the other cases. 



(<7theo) 




c 


n 


average error 


<?mo) 


0.25 






0.08 


(«<o) 




0.34 




0.08 


(q»v) 


0.18 




0.07 


0.07 


(QCv) 




0.26 


0.09 


0.03 



(H) and (I) since the outflow in the case (G) (n e = 4) 
and in the case (H) and (I) are similarly small as we have 
discussed in the Sec. IVII A[ This assumption is also sup- 
ported by the fact that pedestrians stands between two 
former pedestrians to see their way clear in the evacua- 
tion situation, which is achieved by the Triangle-lattice 
neighborhood (n e = 4). 

Table HI] shows the parameters calculated by the 
method of least squares and the average errors for each 
outflow. The errors of (q^o), (q(o), an d (<?c?) are iar g er 
than 0.05. This means that only considering one factor, 
i.e., "the difference of the strength of conflicts according 
to the number of involved pedestrians" or "the decrease 
in walking speed by turning" , is not sufficient to repro- 
duce the realistic outflow. By contrast, the error becomes 
smaller than 0.05, when both the frictional and turning 
functions are introduced. This result verifies that the 
introduction of the frictional function and the turning 
function enables us to obtain the more realistic pedes- 
trian outflows, theoretically. 

The average pedestrian outflows of the theoretical 
analysis are compared with those of the experiments in 
Fig [T7] We can clearly see the four remarkable points 
discussed in Sec. IVII Al (q( v ), which includes the effect 
of both the frictional function and the turning function, 
agrees with the experimental results well except (B). On 
the contrary, (q^o) is not correspond to the experimental 
results well. The effect of turning is not applied to (g^o), 
thus its values are same if n e is same. Moreover, it does 
not use the frictional function either, so that it is impos- 
sible to change its value for n e appropriately. Therefore, 
we confirm that the frictional function and the turning 
function have an essential role to calculate the realistic 
pedestrian outflow again. 

The outflow of the experiment (B) is extremely large 
and the result of the theoretical calculation does not 
agree with it even if the two functions are introduced 
because of the zipper effect. Since the pedestrian intel- 
ligence enables them to achieve the smooth mergin g, i t 
should be modeled to reproduce the zipper effect [34| . 
However, when we consider the normal evacuation and 
the lined evacuation whose lines are more than two, the 
effect of the pedestrian intelligence is possible to neglect, 
and also the evacuation in two lines (B) is a special case; 
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FIG. 17: The pedestrian outflows from the experiments and 
the theoretical calculation. The parameters in Tab. [TT] are 
used for the theoretical results. We observe that {q( v } agree 
with the experimental results well except (B), while (q M o) does 
not. 



therefore, our model includes sufficient factors for ana- 
lyzing usual evacuations. 



VIII. EXPERIMENT AND ANALYSIS 2 
EVACUATION THROUGH AN EXIT 
WITH AN OBSTACLE 

A. Experiment 

We study the effect of an obstacle put in front of an 
exit in this section. We did the evacuation experiment 
at the NHK TV studio in Japan. Two large walls were 
set up in the studio, and we could adjust the width of 
the exit by moving them (Fig. [15] (a)). We decided it 
as 50 cm. The participants of the experiment were fifty 
women, who were their thirties and forties. We did three 
kinds of experiments as in Fig. 1181 The experiment (a) 
(Fig. [TH] (a)) was the evacuation in a line, which was the 
same kind of experiment as the experiment (A) in Fig. 1151 
(b)-(A). The experiment (b) (Fig. [H](b)) was the normal 
evacuation, whose initial condition was the same as the 
experiment (H) in Fig. HH](b)-(H). The initial condition 
and the way of the evacuation in the experiment (c) (Fig. 
[H](c)) was the same as the experiment (b), however, the 
column, whose diameter was 20 cm, was put in front of 
the exit. 

The pedestrian outflows of the three cases are de- 
scribed in Tab. IIIII We calculate them by using (|20[) . 
where i and j are determined as Tab. IIIII by consider- 
ing the same reason as in Sec. IVII Al The outflow in 
the case (a) is the largest, since there was no conflict. 
Comparing the result of the experiment (b) and (c), we 
surprisingly find that the outflow of the experiment (c), 
i.e., the experiment putting an obstacle in front of the 
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(a) One Line (b) Normal (c) Obstacle 

FIG. 18: Schematic views of the experiments, (a) "Going 
through the exit in a line". (Same kind of experiment as 
Fig. [15] (b)-(A).) (b) "Normal evacuation". (Same kind of 
experiment as Fig. [15] (b)-(H).) (c) "Putting an obstacle in 
front of the exit". The width of the exit is 50 cm and the 
number of the pedestrians is fifty. We consider that the four 
pedestrians try to move to the exit at the same time in (b), 
while three pedestrians do in (c) since there is an obstacle. 



exit, is larger than that of the experiment (b), which is a 
normal evacuation. The P value of the difference of the 
average pedestrian outflow between (b) and (c) is calcu- 
lated as p < 0.01, which verifies that the difference is 
significant. 



B. Theoretical Analysis 

We explain this phenomenon by our theoretical cal- 
culation. Watching the video of the experiment, we see 
that the obstacle blocks the pedestrians moving to the 
exit. In the Sec. IVII Bl we have verified that there are 
approximately four pedestrians move to the exit at the 
same time, i.e., n e = 4, in the normal evacuation (Fig. 1181 
(b)). We assume that an obstacle blocks the pedestrians 
moving, so that the number of the pedestrians moving to 
the exit at the same time decreases by one. Therefore, 
we have modeled the obstacle by blocking one neighbor- 
ing cell as in Fig. QjJ] (a) and consider that n e = 3 in the 
experiment (c) (Fig. [TH](c)). This modeling is justified 
by the experimental result that 35 % pedestrians went 
through the left side of the obstacle, whereas 65 % of 



TABLE III: The pedestrian outflows obtained from the ex- 
periments. The data show that the outflow increases if we 
put an obstacle in front of the exit in a proper way. We verify 
that there is a significant difference between the case (b) and 
(c) by calculating P value, which is obtained as p < 0.01. 



Case 


n e 


6»e 


<<7e X p) 


N 
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3 


(a) One line 


1 


(0°) 


3.23 


3 


4 


50 


(b) Normal 


4 


(90°, 30°, 30°, 90°) 


2.80 


6 


1 


47 


(c) Obstacle 


3 


(90°, 30°, 90°) 


2.92 


6 


1 


47 



TABLE IV: Parameters and the average errors for the exper- 
iment (a) to (c) from the method of least squares, /_t, £, and 
r\ calculated by the least square method are described in the 
second, third, and forth columns, respectively. In the fifth 
column, the average error is described. We observe that there 
are errors order of 10 -2 for (g M o), (gco)i an d while (?<;,,) 

agrees with the experimental data very well. 







c 


n 


average error 




0.23 






0.05 


(<7C0> 




0.27 




0.04 




0.23 







0.05 


(907 ) 




0.22 


0.09 






them went through the right side of it, that is, the ratio 
between the pedestrians going through the left side and 
right side of the obstacle is approximately 1:2, which 
corresponds to the number of the left and right neighbor- 
ing cells in regard to the obstacle cell in Fig. [19] (a). 

Now, we obtain the parameter (3 from the result of the 
experiment (a) as 



1, 



0°)) 



3.23 x Alx At, 
0.97. 



(24) 



Then, we calculate [i or £, 77, and the average errors (|21|1 
by the method of the least squares. 

The results are shown in Tab. IIVI We see that the 
average pedestrian outflows calculated by using the fric- 
tional function <f>^ and the turning function t, i.e., (q^ri), 
reproduces the experimental outflow very well since there 
is no error. However, there are average errors whose or- 
ders are 10 -2 for other three theoretical outflows, i.e., 
(<?mo)j (<?Co)j an d (ifio) ■ The frictional parameter does not 
change the strength of the clogging against n e ; therefore, 
the difference of the outflow in the case (b) and (c) is not 
appropriately reproduced. When we adopt the frictional 
parameter, we also find that the turning function can- 
not work effectively since 77 = in the case (c^i,) (Tab. 
lIVp . On the contrary, the frictional function can adjust 
the strength of clogging adequately, thus, the result of 
the experiment is obtained if it is used with the turn- 
ing function. This result also verify the importance of 
the frictional function and justifies our assumption that 
the obstacle increases the pedestrian outflow since it de- 
creases n e , which is the number of pedestrians moving to 
the exit at the same time. The turning function plays an 
important role to represent the position of the obstacle. 
It enables us to distinguish which cell is blocked by the 
obstacle, the cell whose incident angle is 30° or 90°. If we 
do not use the turning function, the effect of the obsta- 
cle's position, which is studied in the next section, is not 
represented properly, so that the average error becomes 
large in (q co ). 

Figure [20l shows the experimental and theoretical av- 
erage pedestrian outflows. The well correspondence be- 
tween the experimental flows and {q^ v } is observed. How- 
ever, (<7 M o) is almost same for the case (b) and (c), and 
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(a) Shifted-Obstacle (b) Center-Obstacle 

FIG. 19: Schematic view of the exit cell and its neighbors 
when an obstacle is set up. (a) The obstacle is shifted from 
the center, (b) The obstacle is at the center. We consider the 
incoming pedestrian flow from the diagonal neighboring cells 
of the exit cell since the obstacle is much smaller than the one 
cell. 
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does not agree with the experimental results. 



FIG. 20: The pedestrian outflows from the experiments and 
the theoretical calculation. The parameters in Tab. [TT| are 
used for the theoretical results. We observe that {q^ v } agree 
with the experimental results well, while (g M o) does not. 



C. Effect of the Position of the Obstacle 

We consider how the position of the obstacle affects on 
the pedestrian outflow in this section. Fig. [19] (b) repre- 
sents the case that the obstacle used in the experiment 
(c) in Fig. [TH] (c) is set at the center of the exit. We 
assume that the pedestrians also move to the exit cell 
from the two diagonal neighboring cells since the diame- 
ter of the obstacle (20cm) is much smaller than the size 
of the cell (50 cm) . We calculate the pedestrian outflow 
in this case by using (q^) with the parameters in Tab. 
HVl which are (3 = 0.97, £ = 0.22, and 77 = 0.09 as follows: 

(tenter) = (q(n e = 4, e = (90°, 45°, 45°, 90°))) ^ 



This value is not only smaller than the outflow of the 
experiment (c) , which is a evacuation with an shifted ob- 
stacle, but also the outflow of the experiment (b) , which 
is a normal evacuation. The result indicates that the ob- 
stacle does not always increase the pedestrian outflow. If 
the position of the obstacle is not adequate, it does not 
block one pedestrian moving to the exit cell; therefore, 
the outflow is not improved. Moreover, if the obstacle 
is at the center of the exit, pedestrians whose incident 
angles are 0°, i.e., who are going to go through the exit 
straightly, have to detour the obstacle, so that their walk- 
ing speeds decrease and the outflow decreases. Thus, the 
obstacle should be set up the place where it blocks the 
pedestrians interrupting from the side and does not block 
the pedestrian walking straight to the exit. We have dis- 
covered this phenomenon since we consider the effect of 
conflicts and turning by the frictional function and the 
turning function. 



IX. CONCLUSION 

In this paper we have introduced the frictional func- 
tion and the turning function to the floor field model. 
The frictional function is the improvement of the fric- 
tion parameter and changes its value against the num- 
ber of pedestrians involved in a conflict. The turning 
function represents the decrease in walking speed when a 
pedestrian change their orientation. Since conflicts and 
turning are the two major phenomena observed at an 
exit in evacuation situation, introduction of the two func- 
tions enables us to calculate realistic pedestrian outflows 
through an exit. 

We have done the two evacuation experiments. One 
is the evacuation in lines in various patterns, and four 
interesting phenomena have observed. First, the pedes- 
trian outflow decreases when the number of pedestrians 
moving to the exit increases since the number and the 
strength of conflicts increases. Second, it also decreases 
when pedestrians need to change their orientation at the 
exit. Third, when pedestrians evacuate in two parallel 
lines, they are intelligent enough to avoid conflicts and 
smoothly merge into one line. Forth, we have discovered 
that approximate four pedestrians are trying to move to 
the exit at the same time in the normal evacuation when 
the width of the exit is 50 cm. 

By using the frictional function and the turning func- 
tion, we have succeeded to reproduce the characteristics 
described in above except the third one and obtained the 
more realistic figure of the pedestrian outflow through an 
exit, which corresponds to the result of the experiments 
very well. The pedestrian outflow using the friction pa- 
rameter or excluding the effect of turning does not agree 
with the experimental results well; therefore, the fric- 
tional function and the turning function are necessary 
for realistic pedestrian outflows. 
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The other experiment is the evacuation through an exit 
with an obstacle. We have clearly shown that the pedes- 
trian outflow increases by putting an obstacle in front 
of the exit from our experiments. Our assumption that 
the pedestrian outflow increases since the obstacle de- 
creases the conflicts at the exit by blocking the pedestri- 
ans' movement is verified by theoretical calculation using 
frictional function. We have also discovered that the out- 
flow depends on the position of the obstacle. When the 
obstacle is shifted from the center of the exit, it blocks 
the interrupting pedestrians form the side and the out- 
flow increases. On the contrary, when the obstacle is 
put at the center, it blocks the pedestrians trying to go 
through the exit straightly and the outflow decreases. We 
have succeeded to obtain this result by introducing the 
turning function. 

While the real pedestrians are intelligent enough to 
merge smoothly in the evacuation in two lines as we dis- 



covered from our experiments, pedestrians in our model 
are not intelligent enough to do so. Thus, developing a 
model which includes such pedestrian intelligence is our 
future work. The effect of a size of an obstacle and the 
width of the exit should be also studied in detail by using 
both the frictional function and the turning function in 
the near future. 
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